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A new approach called the Schrodinger Collective Intrinsic Model (SCIM) has been developed to 
achieve a microscopic description of the coupling between collective and intrinsic excitations. The 
derivation of the SCIM proceeds in two steps. The first step is based on a generalization of the sym- 
metric moment expansion of the equations derived in the framework of the Generator Coordinate 
Method (GCM), when both Hartree-Fock-Bogoliubov (HFB) states and two-quasi-particle excita- 
tions are taken into account as basis states. The second step consists in reducing the generalized 
Hill and Wheeler equation to a simpler form to extract a Schrodinger-like equation. The validity of 
the approach is discussed by means of results obtained for the overlap kernel between HFB states 
and two-quasi-particle excitations at different deformations. 

I. INTRODUCTION 

The Generator Coordinate Method (GCM) is a very useful approach to study large-amplitude collective modes. A 
review on the GCM and quantized collective motion is given in reference [lj . It is widely used for the study of ground 
state properties, and low-lying collective states in even-even nuclei. For instance, in Ref. 0| the Hill- Wheeler equation 
is solved for the axial quadrupole collective coordinate, where the basis states are angular-momentum and particle- 
number projected Hartree-Fock-Bogoliubov states using the Skyrme energy density functional. Improvements have 
been made recently to generalize this study by taking into account the triaxial degree of freedom, with first results 
reported for 24 Mg [3j. State of the art configuration mixing calculations based on the GCM, including variation after 
projection effects and projection of angular momentum and particle number for the effective interaction D1S have 
also been performed by the Madrid group [1, H[. In Ref. Q, shell closures have been studied in the N — 32,34 
mass regions and in [5j shape coexistence near neutron number N = 20 has been analyzed. When used with the 
Gaussian Overlap Approximation (GOA) and employed with the full quadrupole coordinate the GCM leads to a five 
dimensional collective Hamiltonian. This model has been extensively put to the test recently in Ref. [(| @] , where 
predictions for yrast states up to 6 + , and non-yrast 0^, 2% , 2% states have been compared with experimental data. A 
time-dependent version of the GCM has also been used to study the fission process [8, 9]. In these latter works, the 
collective dynamics is derived from a time-dependent quantum-mechanical formalism where the wave function of the 
system is of GCM form, and for which a reduction of the GCM equation to a Schrodinger equation has been made 
by means of the usual techniques based on the Gaussian overlap approximation. 

In all the aforementioned works the adiabatic assumption has been made and no intrinsic excitations have been 
taken into account. Collective vibrations have been considered to be completely decoupled from intrinsic excitations. 
In the seventies, a few pioneering works were published, in which the usual GCM was extended to include 2 quasi- 
particle (qp) excitations (GCM+QP) 0, El- In Ref. [ll|, the method was applied to A ~ 50 nuclei, for which 
excitation spectra display both collective and non collective features. In fact, calculations showed that vibrational 
degrees of freedom may be important to describe the 2 + and 4 + states, while 6 + states may be dominated by 2-qp 
components. In these calculations, the whole pf shell was used as a basis assuming 40 Ca as an inert core, and a 
modified version of the Kuo matrix elements was considered for the residual interaction. A GCM+QP approach is 
well suited for the study of many phenomena, for which the coupling between collective and intrinsic states may play 
a role: i) shape coexistence where K = individual states could be coupled to low-lying vibrational collective states, 
ii) backbending phenomena in rotational bands, iii) decay of superdeformed states to normally deformed states and 
iv) non adiabatic effects during the fission process. 

Where the fully microscopic treatment of the whole fission process from the initial stage of the fissioning system up 
to scission is concerned, GCM based approaches are the best suited. They enable one to treat the fission dynamics as a 
time-dependent evolution in a collective space. Adiabatic calculations have shown that the dynamics is very important 
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for fission fragment distributions Q: it is responsible for the large widths of the distributions. These results have 
validated the adiabatic hypothesis made as a first approximation for the description of low-energy fission. In fact, 
nuclear superfluidity induced by pairing correlations, in addition to strongly influencing the magnitude of the collective 
inertia, is responsible for a strong rearrangement of the individual orbits with deformation. However, there are some 
experimental indications that pairs are broken during the fission process. For instance, manifestations of proton pair 
breaking are observed in 238 U and 239 U nuclei for an excitation energy of 2.3 MeV above the barrier: first the proton 
odd-even effect observed in the fragment mass distributions decreases exponentially for an excitation energy slightly 
higher than 2.3 MeV [l2| and then the total kinetic energy drops suddenly (l3l . fl4| . Some theoretical calculations 
have studied the dissipation during the fission process, most of them are based on a semi-classical formalism such as 

, or Hamiltonian equations with one body dissipation and two-body viscosity fl8l . [l9| , 
|. Here we aim at developing a microscopic non-adiabatic Schrodinger equation, in 
order to obtain a microscopic description of the coupling between collective and intrinsic excitations. We use in the 
following the abbreviation SCIM which stands for Schrodinger Collective Intrinsic Model. The newly developed SCIM 
formalism is presented in Section [H] In Section ITTT| the calculation of the N-dimensional overlap kernel between and 
2qp HFB states is shown, and the dependence of its moments on the deformation is analyzed in details. Section IIVI 
is then devoted to the derivation of the inverse of the overlap kernel. The Schrodinger equation is derived in Section 
Ivl and conclusions are drawn in Section I VII 



Focker-Plank equations ll 
or Langevin equations [2 




II. DESCRIPTION OF THE FORMALISM USED TO CONSTRUCT THE SCIM 



Our derivation of the SCIM proceeds in two steps as described in the two following paragraphs. The first step 
is based on a straightforward generalization of the symmetric moment expansion of the equations derived in the 
framework of the GCM [22]. We refer to the extensive review of the GCM given in [23j . where most of the material 
and references related to this subject are available. The second step consists in reducing the equation obtained in 
section III Al to a simpler form which will be the starting point to extract a Schrodinger- like equation. 



A. Generalization of the Symmetric Moment Expansion 



In order to take into account the coupling between collective and intrinsic degrees of freedom we use the ansatz 

I*} = / dqfo(q)\Mq)) + E / dqfMlMl))- (1) 

In the present study, &o(q) is the ground state at the q value of the collective variable and $i(<7), with i ^ 0, are 
the associated 2-qp excitations whose precise definition and properties are given in section 1111 Al Following the usual 
procedure, the weight functions fi(q) are found by requiring that the total energy (^>\H\^>) / (^l 1 ^) calculated with 
the function defined in Eq. Jl} be stationary with respect to arbitrary variation 5f*, which leads to the well known 
Hill- Wheeler equation 

E / # ( H v (?« «0 - EN * («' <?')) fi = ( 2 ) 

3 J 

where the Hamiltonian and overlap kernels are defined in terms of matrix elements such as 

Nuhd) = (Mq')\^(q)). 

Eq. @ could be the startin g po int to derive a local expansion by means of a Taylor expansion of the fi(q') around 
q, as performed in reference [24]. However we prefer here the symmetric expansion [22j which, in our point of view, 
provides a natural way to expand on the non-locality. Its derivation is very simple if we express the variational 
principle after performing the change of the variables q — {q + q')/2 and s = q — q' in the expression of the total 
energy. We only give here the resulting Hill- Wheeler equation: 
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dse- p / 2 [H{q+ -,q- -) - EN(q+-,q- -)] e isP / 2 f(q) = (3) 
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with 

Let us mention that in all the equations of this formalism, we set h = 1. At this stage it is important to mention that 
Eq. is obtained after performing successive integration by parts and consequently the form given here supposes 
that the corresponding contributions of the integrated terms (surface terms) vanish at the boundary of the integration 
domain. Note that this is the case in all applications of the GCM in spectroscopy studies [(| 0], as well as in nuclear 
fission studies @, In the following, we assume that we are in such a situation and pursue our derivation with Eq. 
© as given above. Now, the series expansion of the exponentials is inserted in Eq. ^ and terms of the same order 
in s are collected together. After integration with respect to s one is led to a symmetric moment expansion that 
includes the coupling between collective and intrinsic variables. Following the notations in reference [23j, except for 
the imaginary number 'V, the moments of any operator A is defined as: 



d SS "A(q + -,q--) (5) 

and symmetric ordered products of operators A^ n ' (q) and P as 

[AW(q)P] (n) = ±Y, C * pn ~ kA(n) (d pk - ( 6 ) 
1 k 

With the procedure described above and the notations just given, Eq. ([3]) is transformed into 

£ I ( [ H W (q)P] {n) -E [NV> [q)P] in) ) m = 0- (7) 
n 

Eq. is the compact form of a set of coupled equations for the different components fi(q). Notice that the moments 
occurring in its definition are square matrices whose dimension depends on the number of excitations introduced in the 
description. It is worth pointing out here that another difference with the usual approach (no qp excitations) is that 
all moments, odd or even, must be included in the summation. In order to give some more information about Eq. ([7]) 
we refer to Section Hill where we elaborate on the properties of the selected GCM collective excitations in the present 
work. There, it is shown that, with such a selection, the even moments and odd moments of Hermitian operators are 
represented by real symmetric matrices and imaginary anti-symmetric matrices respectively. As a consequence the 
operators in Eq. {7} are Hermitian. Finally, since the Hamiltonian is time-reversal invariant, it is easy to check that 
this operator is also invariant if collective and intrinsic coordinates are time reversed as well. We will comment later 
on the effect this result has on the Schrodinger equation we derive in Section fVl 



B. Reduction of the Symmetric Moment Expansion to a Schrodinger like equation 

Eq. ([7]) is an infinite series which is an exact expansion of the Hill- Wheeler equations in terms of local operators. 
It can be transformed into a local Schrodinger equation by inverting the expansion of the overlap kernel. The latter 
can be written formally as 

i%)=£^"%)P] (n) = (N^(q)) + N^(q). (8) 

n 

In the following we use the definitions 



N^(q) = J^ 2 (q)^/N(0)(q) 
(N^(q)) + = y/W°m(jV 2 (q)) + , (9) 



with 

J(q) = I + u(q) 



u(q) - 1 r^[jyW( t -) f |W_ J (io) 
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These notations require further explanation. We use the property that any Hermitian operator A can be written in 
the form A = S + S, with an operator S which is not necessarily Hermitian. Then, for the sake of simplicity, the 
notation S = A 1 / 2 is used although it is not strictly speaking the square root of A unless S + = S. 
We proceed further by defining a new set of collective wave functions, {gi(q)}, through the matrix relation, 

/(g) - N-^iqMq) = {N^m-^J-^^giq) (11) 

to arrive finally from Eq. {7} at the local Schrodinger equation 

( {j~i/2{q)) + \Y , 1 = [ff (n) (g)Pl (n) - 1 lj-i/a(g)--E|g(g) = 0, (12) 
where we have defined the operator 

J_ 1/a («) = J- 1/2 (q). (13) 

Eq. is formally the same as the one given in [24| but differs from it in many respects. The definitions of 

the moments and the operator u(q) are not at all the same. Furthermore, the operators occurring in its definition 
are represented by matrices including, as already mentioned, those corresponding to odd moments. Notice that the 
normalization 

' dq(g(q)) + g(q) = l (14) 

guarantees the normalization of the wave function defined in Eq. ([T]). In this approach, the main difficulty is to 
determine the operator J_i/2(<z), which, according to the definitions given in Eq. (|10[) . is solution of the equation 

{■J- 1/2 (q)) + (I + u(q))J_ 1/2 (q) = I. (15) 

In order to reduce its complexity, one is led to assume that the series expansion representing the overlap kernel 
converges rapidly and can be truncated after the second-order moment. That is to say, we approach Eq. (TPoj) by the 
simpler one 

v-^ + [ I+ jmtfi" < *'M m + \[^w<^ m ]j- 1 Mt) - 1. del 

In section IIVI we study Eq. (|16[) in some detail and explain how it can be solved with parametrizations of the form 

4 

71=0 

The quantities jr n \ (q) are unknown matrices which are determined by solving Eq. (I16p . 

Once the inverse J_ 1/2(9") is known it can be inserted in Eq. (|12[) . and by means of the formulas given in the 
Appendix [^1 there is no difficulty in finding a general expansion in terms of symmetric ordered products of operators. 
In other words one can express Eq. (fT2)) in the form 

J2([S {n) (q)P] in) -E)g(q) = 0. (18) 

n 

Then a second order differential Schrodinger equation is derived by limiting the summation to n = 2. However, for 
obvious practical reasons, one does not extract the S^ n \q) with the full expansion of the Hamiltonian kernel but with 
a truncation similar to the one introduced in the case of the overlap kernel. Thus, in all that follows, our derivation 
of a Schrodinger equation relies on the approximated expression 



( J-i/2(q)T 



(H^(q)+[H^(q)P] 



(1) 



l[ff( 2 )(»)pi (2) ) . 1 



j-i/2(q)-E)g(q)=Q, 



(19) 
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with .7-1/2(9) satisfying Eq. (jTHJ) . 

Such an approach assumes that the series expansions of the overlap and Hamiltonian kernels converge rapidly. If the 
moments were independent of q or if we could neglect their variations, Eq. ([7]) would simply be a power expansion in 
the collective momentum P. Then one could conclude naturally that this approach is limited to the study of collective 
motions with momenta that are not too large, or equivalently at moderate energies. Notice, that if we approximate 

the kernels with a Gaussian [23| of width a(q) the expansion (Eq. ([7])) becomes a series with terms [a n {q)P] [n) . In 
that form it appears that a rapid convergence of Eq. ([7]) requires that the collective wave functions, as well as the 
moments, are slowly varying functions in a range given by the width of the Gaussian. These arguments are only 
qualitative in the sense that they do not give precise information regarding the energies up to which one can expect 
the approximation to be reasonable. In relation to this discussion we refer to the work by Bonche et al. [HJ which 
contains a detailed comparison between two calculations of large amplitude quadrupole dynamics in 194 Hg. One is 
based on a direct resolution of the Hill- Wheeler equation (HWE) as extracted with the GCM, while the other uses a 
collective Schrodinger equation (CSE) resulting from a local expansion of the HWE into a second order differential 
equation. Although their approach to the CSE approximation is not rigorously identical to the one described here, 
we think that one can extract general deductions from their work. We will only quote some of their conclusions 
here. The GCM and the CSE reproduced the energies of most collective states fairly well. However one observes 
significant differences in the collective wave functions at excitation energies around 6 MeV and above. More precisely 
the corresponding HWE collective wave functions display more and more irregular structures and rapid variations 
that a CSE approach is not able to reproduce. We do not deny the fact that it is important to have a fair description 
of the collective wave functions, but we do make the observation that for excitations above 4 or 5 MeV one also has 
to worry about their possible coupling with qp excitations, which is the main motivation of this work. 

At this stage one needs some more information on the kernel overlap to proceed further in the reduction of the 
equations given above. It is the purpose of the next section to present a numerical study of this kernel in a wide range 
of deformations that would be encountered in a microscopic description of the symmetric fission of 236 U. 

III. THE N-DIMENSIONAL OVERLAP KERNEL 

A. Selection of the 2-qp generating wave function 

In this paragraph, a set of generating wave functions is constructed by means of the constrained Hartree-Fock- 
Bogoliubov method as described in [8, 9]. We limit this study to the case of one generator coordinate q associated to 
the total quadrupole deformation of the even-even nucleus 236 U. Furthermore, the calculations are restricted to axial 
symmetry. As a consequence the qp are characterized by their projection K of the total angular momentum on the 
symmetry axis. We denote the ground state (or Bogoliubov vacuum) at deformation q by |<&o(<z)) an d recall here its 
definition in terms of the qp destruction operators r)i(q) and their time reversal conjugate Tji(q) = T + r]i{q)T, where T 
is the time reversal operator: 

l*o(?))=n»?*(9)'B(9)l )- ( 2 °) 

i>0 

Furthermore we restrict our calculations by introducing another symmetry which is associated with the invariance 
of the Bogoliubov fields when performing a reflexion with respect to the xQz plane. In that case the parity tt is an 
additional quantum number and the qp operators are then characterized by k = {n^, 7Tfc, K^} where rik is an index 
which distinguishes single qp states inside the same block {irk,Kk}. At each deformation we can then define a set of 
2-qp excitations 

v+MK(i)\Mq))- (21) 

This study is limited to the description of the coupling of the ground state with 2-qp excitations. Consequently, 
since the Hill- Wheeler equation conserves the total projection of the angular momentum on the symmetry axis (z) 
and since the ground state is a state K — 0, it is clear that the quantum number of 2-qp excited states must satisfy 
Ki t + Ki 2 = 0. In addition, by using similar arguments with the time reversal transformation, we arrive at the 
conclusion that the qp excitations must be symmetric under time reversal. Thus the TV selected configurations at 
deformation q are written in the form: 

' 1*0(9)) 

< |*i(9)) = ^[l + *<x«(^-l)]Kt(«)<(?) (22) 

I $0 (<?)>, i = l,N-l. 
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Notice that for the 2-qp states, the index i denotes the couple (£1,^2). For simplicity of notation we have not added 
another index to indicate whether the excited configurations are built with two neutrons or two protons qp. Finally, 
for sake of convenience and since the 2-qps are m t ne same block {iTi,Ki} = {^iu 2 ) ^h(2)} ' we a ^ so use m 

the following the notation i = {K?*, (ai)} where the letter which distinguishes the couple (71^,71^) in the same 
block {iTi^Ki} is used if necessary. These wave functions which are not eigenfunctions of the proton and neutron 
number operators N p , N n should in principle be projected on the subspace of states with good particle numbers. The 
expression of the ove rlap kernel is readily expressed with the usual techniques to calculate the overlap between two 
Bogoliubov vacuum |27j. The difficulty in projecting is essentially numerical since the projection method requires 
repeating the calculation of the overlap at different angles in order to calculate an integral with a discretization method 
[28j. Such calculation, if feasible in the study proposed here, would represent a considerable amount of computational 
work and numerical studies. At this stage it is worth pointing out that a description of the fission process from the 
first barrier to scission involves a wide range of deformations typically in the interval [50 b, 550 b]. Furthermore, an 
accurate calculation of the potential energy surface, in particular close to the scission point, necessitates the use of a 
very large dimensional space of harmonic oscillator basis states. As an example, we indicate that in a one-center basis 
the size of the if = 1/2 block is [150 x 150]. In view of these remarks, we find some justification in presenting this 
formalism with unprojected wave functions. We stress that until now our formalism does not depend on the particular 
model used to calculate the kernels. Finally, let us conclude this section with a different issue which is still related to 
this question of particle number. It concerns the restoration of the correct average values of the number of protons 
and neutrons of the Hill- Wheeler solutions. This problem is studied in details in reference [25| where the authors 
introduce in the Hill- Wheeler equation the kernels calculated with the two operators N p , N n and their associated 
Lagrange parameters. The implementation of this method in our formalism is straightforward. 



B. Calculation of the overlap kernel 

1 . Generality 

Calculations presented in this paper correspond to different excitations along the symmetric fission barrier in 236 U. 
The HFB equations are solved by expanding the single-particle states onto an axial harmonic oscillator (HO) basis. 
A one-center basis with N = 14 major shells has been used. The calculations are performed with the Hartree-Fock 
+ BCS method and a constraint on the quadrupole operator. They have been restricted to axial and left-right 
symmetries. The full Bogoliubov approach has not been considered here for reasons of computational time only. 
Because of the symmetries imposed in this calculation, the Bogoliubov or (HF+BCS) transformation can be chosen 
to be real. As a result all quantities occurring in the formalism are real. This is the case in particular of the overlap 
matrix which is denoted by 

N ij ( q ,q') = {^ i (q)^ j ( q ')}. (23) 

Calculations of this matrix have been performed at different elongations along the symmetric barrier. The formalism 
used to extract these matrix elements is given in detail in Appendix [Bj In all what follows the results are presented 
in the variables q and s defined in the first paragraph of Section |TTJ More precisely at each q given in barns, we 
plot the matrix elements Nij(q + s/2,q — s/2) as function of s running in intervals wide enough that the overlaps 
Nij(q+s/2, q— s/2) drop to zero. The {q} have been chosen in the range [20 b, 350 b] with a step of 10 b up to q = 100 b 
and a step of 25 b beyond. Let us mention that deformations close to 160 b, namely q = 150 b and q = 175 b, are not 
discussed here, since they correspond to a junction between two valleys having different hexadecapole deformations. 
Such a study requires to treat explicitly the hexadecapole degree of freedom [23 . |3C| . 

The symmetric barrier is plotted in Fig. [TJ The ground state is located near 30 b, the first barrier at 50 b, the 
second well at 80 b and the second saddle point at 160 b. 

The intrinsic levels chosen to construct the 2-qp excitations are the lowest-energy proton = 
1/2 - , l/2 + , 3/2~, 3/2 + , ... qp levels. Among them, the selected intrinsic excitations are those which minimize the 
deviation from the average proton and neutron numbers. Results are presented in the present work only for proton 
excitations and similar results and conclusions are expected for neutron excitations. 

Some comments are in order concerning phase problems that we encountered in the construction of these 2-qp 
excitations. In fact the qp states or operators r\i which result from the diagonalization of the Bogoliubov Hamiltonian 
are defined up to an arbitrary phase which in the present case is real, i.e. ±1. The phases of two conjugate qps 
rp[,r]i are then the same by construction and consequently the Bogoliubov vacuum defined above does not depend 
on these phases. The same conclusion applies to the case of two conjugate qp excitations \$i(q)) = f^?^ !$(<?)) 
but it is not the case for all other 2qp excitations. It is worth stressing that if we were using the qp as they come 



7 



-1760 



> -1765 




-1780, 



50 100 150 200 250 300 350 400 
q(b) 

Figure 1: HFB potential energy along the symmetric fission barrier in 236 U 



out of the diagonalization of the Bogoliubov Hamiltonian, the corresponding non-diagonal matrix elements of the 
kernels can become discontinuous at random as function of the deformation, which would prevent the extension of 
the GCM with those excitations. In order to solve this problem, the single-qp excitations used in the construction of 
the 2-qp excitations have been followed continuously as a function of deformation. This is achieved with the following 
simple procedure. Any single-qp state \i,q + Sq) at deformation q + Sq is associated with \i, q) by maximizing the 
single-particle overlap \(i,q\k,q + Sq)\ among all k in the 1-qp spectrum at deformation q + Sq. This overlap is 
nothing but the anti-commutator {77+ (q), r/k(q + Sq)}. Denoting by k the qp which achieves this maximization, we set 
\i, q + Sq) = (fii(q)\k, q + Sq) where a phase given by 



(i,q\k,q + 5q) 
\(i,q\k,q + Sq)\ 



(24) 



has been introduced. It is clear that with such procedure the matrix elements of the overlap kernel become continuous 
functions of the deformation. Let us mention that the problem of the sign of the overlap has been recently addressed 
by L. Robledo in [32] for general HFB wave functions that do not have any spatial symmetry (triaxial) and also 
breaking of the time reversal symmetry. 

The evolution of the energy of 2-qp excitations is plotted in Fig. [2] as a function of the deformation around q = 70 b. 
It shows that the lowest excitation is located around 3 MeV. No significant difference is observed between proton and 
neutron excitations. 

Let us mention that we need not limit ourselves to excitations that conserve the average number of particles, but as 
a result the number of excitations increases rapidly when we allow the mean particle numbers JYj and Z to vary by 1 
particle around their average value. For instance for an energy below only 5 MeV and with a mean particle number 
in the range 92 — 1 < Zi < 92 + 1, 144 — 1 < Ni < 144 + 1, we already find 6.5 proton 2-qp excitations and 11.5 
neutrons per deformation. 

As we follow the single-qp states by continuity as a function of deformation, levels with the same K and 7r quantum 
numbers can occasionally cross or be pushed back from each other. Such features are common and can occur anywhere 
along the deformation. Such cases are illustrated in Fig. [3] around q = 30 b for K* = 1/2+ . The figure shows that 
a repulsion occurs between the two first levels at a deformation q r — 29.5 b, defined as the deformation where the 
difference in energy between the 2 levels is minimal, whereas crossings are observed at q = 34 b and q = 40 b. Note 
that, unlike crossings, repulsions correspond to a mixing of the two incoming levels. 



C. Numerical results for the matrices Nij(q + s/2, q — s/2) of the overlap kernel 

As already mentioned in (25| . an accurate calculation of the complex behavior of the overlap kernel requires a careful 
determination of the HFB states. The overlap kernels are the determinants of the matrices built from the overlaps of 
all pairs of individual orbits: they are very sensitive to the details of those orbits. For this reason one must achieve 
a very accurate convergence of the HFB wave functions if one wants to get reliable values of the overlap kernels. As 
a consequence, all the HFB states here have been obtained with a very high degree of accuracy. More precisely, a 
convergence of 10 -6 fm -3 has been achieved on the generalized density matrix. Let us note that the parameters of 
the Harmonic Oscillators used to develop the HFB states have been optimized in each region at the deformation q 
(i.e. q — q'), and kept constant for the neighboring deformations. For reference we have thought it interesting to 
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Figure 2: Energy of the selected 2-qp excitations around q = 70 b, as a function of deformation. Letters are used to 

differentiate excitations with same . 




show the plot of these matrices as function of s for different values of q even though the quantities which count in 
this formalism are their moments. The latter are presented in section fill C 31 



1. Calculation of the diagonal matrices Nu(q + s/2, q — s/2) 

In Fig. [4]the overlaps Noo(q + s/2, q— s/2) between HFB solutions at different deformations are plotted as functions 
of s for different values of q. 

At first glance, all the overlaps shown in Fig. 3] suggest that they are close to a Gaussian shape, independently of 
the mean quadrupole deformation q. Although the SCIM formalism derived in this paper does not require that the 
overlaps be Gaussian, it is instructive to check the validity of this assumption, as a further motivation for the more 
general moment-based approach we use. The width of the overlap is found to strongly depend on the deformation. 
In general, the width increases with the deformation, with local minima all along the barrier. The extent to which 
the Gaussian shape is a good approximation can be investigated by comparing different ways to associate a width to 
these overlaps. In fact, assuming that they can be represented by a Gaussian 

I(q, q >) = e-^ 2 / 2 ° 2 , (25) 
their moments of even order should verify the relations 

AK%) = V2^a(q) 
N^\q) = (_l)»|^ 2 «( g -)/v-(°%) (26) 
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Figure 4: (color online): Comparison between the overlaps Noo(q + s/2, q — s/2) calculated for different q as 
functions of s for 20b < q < 100b (left panel) and 100b < q < 350b (right panel). 



which provide different ways to extract the width cr(q). The same procedure can be applied a priori to all the diagonal 
terms of N(q + s/2,q — s/2). In Fig. [5] (left panel) a few diagonal elements of the overlap matrix are plotted at 
q = 50 b where a comparison is made between overlaps of HFB vacuum and 2 qp states. The different overlaps are 
very close together, but we notice some difference in the tails. As we will see latter, these small deviations in the 
tails have a significant influence on the value of the second-order moments. Furthermore, we find that level crossings 
do not generate any sizable change in N^\q), whereas repulsions are responsible for a quite rapid but continuous 
change in the overlap. Three different cases are plotted in Fig. \5\ (right panel). First, we observe that the overlap 
(plotted as dashed line) of the excitation built with two if* = 3/2+ qp, one of them being involved in a repulsion 
at q r = 57 b, is not sizably disturbed by the repulsion. Conversely, for the same configuration - one single excitation 
involved in a repulsion at q r = 57 b - the overlap built with Kf = l/2 + qp is sharper (smaller width) than the 
previous one. The most pathological case is found for the overlap built this time with two 5/2 + excitations involved 
in a mutual repulsion at q r — 47.5 b, close to q — 50 b. This overlap even becomes negative in the region close to the 
level repulsion. The comparison with a Gaussian overlap is therefore clearly unjustified here. Let us note that these 
numerical results are in agreement with the analysis of the repulsions we made in Section IHIB II 




Figure 5: (color online): (Left panel) Overlap kernel Nu(q + s/2, q — s/2) for HFB solutions and different 2-qp 
excitations not involved in repulsions (labeled by K*(i)) at q = 60 b, plotted as functions of s. (Right panel) 
Overlap kernel Na(q + s/2, q — s/2) for HFB solutions and different 2-qp excitations involved in repulsions (labeled 
by K v (i), and q r the deformation where the repulsion occurs) at q = 50b, plotted as functions of s. Letters are used 

to differentiate same K v excitations in both figures. 

The widths, obtained from Eq. (|26|) using the moments of zero, second and fourth orders are plotted in Fig. [6] as 
functions of the deformation for the overlaps N 00 (q + s/2, q — s/2) and Nu(q + s/2, q — s/2). 

Particular attention has been paid to the widths extracted from the N^\q) overlaps: calculations of the widths 
have been performed in the range [20 b, 110 b] with a step of 0.5 b. On the right panel of Fig. [5] we observe that 
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Figure 6: (color online) Widths of the overlap extracted from the moments of zero, second and fourth orders, as 
functions of the elongation, in black circle, red cross and blue triangle, respectively. Full lines correspond to the 
widths from Nqq terms and circles to diagonal terms iVi , i ^ 0. Left panel: Results for a large range of 
deformation [0,350b]. Right panel: Widths from Nqq in the region [0,110b]. 



the three widths for Nqq (q) do not differ too much in most cases, which corroborates the fact that the overlaps are 
reasonably close to those from a Gaussian shape. For the lowest deformations, the largest differences occur when the 
width is minimum, at 46 b, 72 b and 102 b with a maximal deviation of 22% at 72 b, between the widths extracted 
from the moments of zero and fourth orders. 

Widths of the overlaps Nu(q), with i ^ are also plotted up to 350b on the left panel of Fig. [51 Since the 
excitations involved in repulsions do not lead to Gaussian shapes, as illustrated in Fig. [SJd), they are not taken into 
account in this analysis. The maximal deviation between the widths obtained from N^(q), N^(q) and N^(q) for 
a given excitation is found to be about 19%. The average of this maximal deviation over all the different excitations 
and over all the deformations is about 5%. However, we emphasize that these deviations will pose no problem in our 
approach since we do not rely on a Gaussian-overlap assumption at all. 



2. Calculation of the non-diagonal terms Nij(q + s/2, q — s/2) 

Non diagonal overlaps Nij(q + s/2, q — s/2) (i ^ j) are plotted in Fig. [71 Noi(q + s/2, q — s/2) are displayed in the 
left panel, and Nij(q + s/2,q— s/2) in the right panel. For the sake of visibility only a few overlaps are drawn. 
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Figure 7: (Color online) Overlaps Noi(q + s/2, q — s/2) (left panel) and Nij (q + s/2, q — s/2) (right panel) at q — 306, 
for different excitations labeled by as a function of s. Letters are used to differentiate same excitations. 



At this point let us emphasize that non-diagonal overlaps have no a-priori reason to be even or odd functions of s. 
However, they must be zero for s = 0, since, at the same deformation, 2-qp excitations are orthogonal to the ground 
state and any other 2-qp excitation. We observe on this plot that in fact they are odd or even functions of s near 
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the origin, depending on the indices i,j of the excitations. Notice that the amplitudes of these non-diagonal overlaps 
are much smaller than the diagonal ones. Furthermore, the overlaps between two excited states are also small in 
comparison to those between the ground state and excited states. The exceptions are found for qp states having the 
same quantum numbers K^(i) — K^fJ), as shown by the blue curve with open circles in Fig. Wp)- This result is not 
surprising since, using the generalized Wick theorem when K n {i) ^ K n (j), the following relation is found: 

xN 0i (q - s/2, q + s/2)N 0j (q + s/2, q - s/2). (27) 
We also mention here that the overlap matrix is not block-diagonal in isospin, since 

Mlp{q,q') = (<S>(qM;vIrt + ^Mq>)) 

is non zero for ^ Tj. 

Conclusions obtained for the diagonal terms of N^°\q) about repulsions and crossings still apply here: non-diagonal 
overlaps are not disturbed by crossing levels but can vary rapidly when a repulsion occurs. In Fig. [5] the overlap 
Nij for a couple of 2-qp excitations i = (21,22) and j — (ji, J2) with = l/2~ is depicted. In this case, i\ and 

12 are both involved in a mutual repulsion at q r — 218 b and — ji- The amplitude of this overlap is much larger 
than that of the standard (without repulsion) non diagonal overlap. Let us mention that this kind of configuration is 
encountered several times along the deformation and corroborates our observation in Section MI Al that the qp states 
i\ and j\ are mixed in the repulsion area. 
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Figure 8: Non diagonal term Nij(q) built with excitations involved in a repulsion at q r — 218b. 




3. Numerical study of the Moments N\ r - (q), n = 0, 1, 2 

According to our approximation (see Eq. (|16j) ). the overlap kernel has to be analyzed through its moments up to 
the second order. At this point it is useful to introduce some important properties of the moments defined in Section 
III Al Their demonstration being straightforward we only list them here: 

- All moments of a Hermitian operator are Hermitian: 

A+ = A=>A^+ = A^\ 

By applying this property to the overlap kernel, whose matrix elements are real for reasons given at the beginning of 
this section, the matrices associated with even-order moments are found to be real symmetric while those associated 
with odd-order moments are imaginary antisymmetric. 

- Finally, let us indicate that the even moments of an operator have the same parity with respect to the time reversal 
symmetry (T) as the operator considered, while the opposite is true for odd moments 

f+Af = eA f+A^f = (-lfeA^, 
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Table I: Matrix elements of N^(q) at q — 60 b. Centers of repulsion q r are mentioned if needed. Letters are used to 

differentiate same excitations. 



where e = ±1. 

With these properties established, we can now analyze numerical calculations of each of the three moments N^j (g) 
for n = 0,1, 2. 

a. The Zero order moment matrix (q) This matrix plays an essential role in the formalism as observed in Eqs. 
(dnj) and (TT9"]) . In fact it is the inverse of its square root which occurs in these expressions and whose determination 
requires in principle the diagonalization of N^°'(q) with a matrix U(q) which depends itself on the deformation q. 
As a consequence, the transformation of the symmetric operators in Eqs. (|16ll and (fl"9|) to the new representation 
necessitates the calculation of the first and second derivatives of U(q) thereby including a number of terms which 
significantly increase the size of the expressions (1161) and (|19p . In light of the numerical calculations presented above, 
we propose instead a set of reasonable approximations to avoid these unnecessary complications. 

First, in Fig. [HI the moment Nqq\(}), as well as its first and second derivatives, all of which intervene in various 
quantities, are shown at different deformations. 
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Figure 9: (left panel) Ng '(q) and its first and second derivatives (right panel) as functions of the deformation. 



The moment generally increases between 8 b and 37 b as a function of the quadrupole deformation. Its first two 
derivatives (g) and (q) are found to be much smaller by at least one order of magnitude. Their values vary 

rapidly in the interval [-0.9,0.5] for N$' (q) and [-0.4 b~\ 0.6 b" 1 } for N$" (q). 

In Table U we give the matrix elements of N^°'(q) for different excitations at q = 60 b. At this deformation, we 
observe that the diagonal elements do not depend too much on the or 2 qp states under consideration: they are 
found to be between 12.0 b to 13.4 b. Such a feature is also observed for most of the deformations. However, let us 
mention that in a few cases, repulsions between single-particle states induce larger variations of the diagonal elements 
(q). This observation suggests we should introduce the diagonal terms N^y(q) defined as 

N^(q) = ±-Tr(N^(<l)) 
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Table II: Matrix elements of normalized first-order moment N^ x \q) for different excitations at q = 60 b. All the 
values are divided by i in order to be real. Centers of repulsion q r are mentioned if needed. Letters are used to 

differentiate same-if^ excitations. 



with N the dimension of the matrix N^°\q). With such a definition, the matrix N^ oS> (q) can then be written in the 
form 

W <%) . + A( S >], A,,", = 

The diagonal and off-diagonal matrix elements of A are all found to be small for all deformations , and therefore the 
inverse of N'°'(q) can be calculated by means of the series expansion of (/ + A(g)) -1 / 2 . As a result of this study, it 
appears that considering this inverse square root as diagonal is a reasonable approximation. 

At this point, we introduce normalized moments for non zero order moments of the overlap kernel since they are 
the quantities of interest in the formalism developed in sections IIVI and [V] 

AK%) = 1 JvW(g) J p > 1. (28) 

Taking advantage of our previous discussion, these moments are calculated in the approximation that the inverse 
square root of N^°\q) is diagonal. Their properties are discussed in the next paragraph. 

b. First order matrix moment N^'(q) The matrix elements of normalized first order moment N^(q) are given 
in Table [TT1 for different excitations at q = 606. Since N^(q) is antisymmetric, the diagonal terms are zero. From Eq. 

|, it is clear that the value of N^\q) depends on the amplitude of the overlap A^- (q). As shown in Table ILT1 the 



absolute value of N^'(q) is generally small although some specific excitations associated with level repulsions lead to 
non negligible values. 

Over the whole range of deformation, the most singular case concerns two excitations involved in a mutual repulsion 
for which \N^\q)\ = 6.9 b. All along the deformation, only 27 matrix elements |iV^(g)| (out of 612) are higher than 
lb. 

c. Second order moment N^ 2 \q) As we will see in Section (jlVj) . the derivation of the inverse of ^1/2(9) = 
(I + ii(q)) 1 / 2 requires the determination of not only the second-order moments, but also its first and second derivatives. 

For the sake of brevity, we only display in Fig. \TUk ) and [TUb) those quantities in the case of the diagonal element 
in the ground state. 

In Fig. [TUk ) we notice that all these quantities, not only vary rapidly with the deformation but they also undergo 
oscillations of very large amplitude. The moment itself varies in between —5 b 2 and —211 b 2 in a range of deformations 
[20 b, 350 b] . As for the first and second derivatives they are contained in the range [—3 b, 6 b] and [—2.1, 1] respectively 
(Fig. HOb)). 

In Table IIII1 diagonal and non-diagonal normalized second-order moments are given for the HFB minimum and 
the selected excitations at q = 80 b. Diagonal terms are negative due to our definition of the second-order moment . 
As expected from our discussion on non-diagonal matrix element of the overlap, most of the non-diagonal terms are 
at least two orders of magnitude smaller than the diagonal ones. It is also seen that they strongly depend on the 

excitations. Very similar features are observed for all deformations up to q = 350 b. However, the maximum value of 

— (2) 2 — (2) 

a non-diagonal term is found to be sup \N>-'(q) \ = 78.4 b z at q = 250 b, which corresponds to 44% of N^ '(q) in this 

region. This case corresponds to i = 12) and j — (ji, J2) excitations for which repulsion occurs between i\ and j\ 
at q r — 242 b and between ii and j% at q r = 265 b. 
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Figure 10: (Color online) The normalized second order moment N^ 2 \q) (left panel) and its first (black filled circles) 
and second (red filled triangles) derivatives (right panel) of the non excited overlap along the deformation. 
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Table III: Matrix elements of normalized second order moment (q) for different excitations at q — 80 b. Centers 
of repulsion q r are mentioned if needed. Letters are used to differentiate same excitations. 



It is worth stressing that perturbations in the overlap matrix due to repulsions, such as the change in amplitude of 
a non-diagonal term as depicted in Fig. [8l or the reduction of the diagonal terms (see Fig. [SJd), disturbs the values of 
the moments, which are the quantities of interest in the present formalism, in a brutal but still continuous way . 



IV. STUDY AND DETERMINATION OF J ±1/2 (q) 

In the previous section we have seen that the moments and their derivatives vary rapidly as function of the 
deformation. As we will see in the following these quantities occur everywhere in the calculation of ^±1/2(9) (with 
J ±1 / 2 (<z) = ^ ±1 ^ 2 (?)), which seriously complicates the derivation of the inverse of the overlap kernel. This question of 
the dependence of the moments on the deformation has been studied in detail in the GOA [l( , and in the symmetric- 
moment expansion for the one-dimensional case of no intrinsic excitations [33l | (although, even in the one-dimensional 
case, our solution below is more complete), but in the context of the present study we have not found any information 
on this subject in the open literature. In the discussion that follows, we will call a symmetric operator of order 

n despite the fact that, strictly speaking, an order cannot be attributed to such an operator since it contains all orders 
up to n. 



A. Discussion and generalities about the inversion of the one-dimensional overlap kernel 



The difficulty mentioned above already occurs in the standard (one-dimensional) approach and consequently the 
discussion is first restricted to this simple case. Let us recall here that when the moments are independent of the 



15 



deformation, the derivation of the operators in question is straightforward. They are given by [23| 

4i /2 (?) = [i + £(?)] ±1/2 ^i±~%) 

with 

W 2N(°)(q) 

This form is consistent with an expansion of the kernel overlap up to second order in the collective momentum P. 
Besides, we are justified in stopping the expansion of j_i/2(s) at the second order in P since, in the usual derivation 
of a Schrodinger equation, one neglects the derivatives of the Hamiltonian kernel of order greater than two. The 
situation is quite different here since symmetric operators of order four give contributions to the expansion of the 
Hamiltonian kernel involving derivatives less or equal than two. In order to go beyond this approximation, successive 
transformations have been performed which allow J(q) to be expressed in a suitable form for our purposes. First the 
normalized moment N@\q) defined above (Section IIII[) is introduced in the definition of u(q) which leads to 

m = _=L=l[iv( 2 )( g -)p] (2i 1 



= A C W (N^ , AT(O)) + I [JV(2) (q) P ] (2) . 

The derivation of this expression is given in Appendix C. The coefficient is a function of the first and second 
derivatives of the moment of zero order. If the latter were constant, then would vanish. After inserting this 
expression in the definition of J(q) given in equation (|10p we obtain 

J(q) = l + a (0) ( g -) + I[AH%)p] (2) 



with 



(q) = --C^(M 2 \N^). 



with 



«(0)x», 2 

According to Fig. [TTJ (q) is not negligible since it can reach 10% at several deformations. It is possible however 
to rewrite J(q) in the following form 

J(q) = y/A (0) mi + «d)(?) + \ [N ( R\q)P] (2 yA {0) (q) 

■a w (q)= -± CW(N(y,A (0) ) 

V A (o)(q) = V 1 + a (o)(q) . (29) 

Coming back to Fig. [TT]we observe that am(q) is at most 4% (at q — 30b) in the whole range of deformation. At 
this stage it seems reasonable to neglect this term but it is clear that, if a better approximation is needed, one can 
iterate this procedure starting with am which gives the quantity Am and so on. If one refers to Eqs. (|16[) and (|19[) . 
one concludes in the light of this study that the entire formalism remains unchanged with a J{q) still given by 

J(q) = 1 + u(q) =l + \ [NW(q)P] (2) , (30) 

provided we substitute in all expressions a renormalized zero moment N^\q) — YV^Iq Au\N^(q) instead of the usual 
moment N^°\q). An index i is used to take into account the number of times the procedure described above has 
been iterated. From now on the notation is kept unchanged, with the understanding that in fact a renormalized zero 
moment is used. 
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Figure 11: Comparison between (N^ 2 \ N^) and first iterated C^ 2 \N^ 2 \ 1 + a) values as a function 

deformation. 



of the 



Thus we are led to determine J±i/2(9) with J(q) given in Eq. (|10p . In fact, we are mainly interested in the 
determination of J_ 1/2(9), since it is this operator we need to express the contribution of the Hamiltonian kernel in 
Eq. (TK)|) . In the one-dimensional case it must be solution of the equation: 



J-l/2(m + \[N^(q)P] i2) }U /2 (q) = 1 
As suggested in reference [23j|, we attempt to solve Eq. (f3"TT) with an ansatz of the form 



■/-172(f) = A_ 1/2 (q) + [B-i/2(q)P] {2) + [C_ 1/2 (q)P 



I (4) 



(31) 



(32) 



After inserting this ansatz in Eq. (|31[) we are led to a complicated coupled set of non linear equations between 
^-1/2(3)1 -B— ill (s) j C- 1/2 (?) an d their derivatives up to high orders. Solving these equations, including the second 
derivative of N^ 2 \q) is not a trivial problem and it would be a considerable amount of work. This is an extensive 
study in its own right that goes far beyond the scope of the present paper. In the following we take advantage of 
the fact that the size of the problem is reduced considerably if we neglect these second derivatives and consequently 
we limit ourselves to finding solutions including only the first derivative of N^(q). Even in this simple case it is 
important to find approximated solutions as a starting point in order to overcome some difficulties in solving Eq. (pH]) . 
A natural choice for these approximated solutions is to approach J± 1/2(9) by the first terms of their series expansion. 
For instance we will define a .7,4,1/2(9) as: 



Ja,1/2{q) = 1 + - a" 2 (?) ~ (1 + «(?)) 



1/2 



where the operator u(q) is defined by Eq. 
operator in the form: 



30]) . By inserting the definition of u(q) in </a, 1/2(9) we readily obtain this 



^,1/2(9) = Ay 2 (q) + [B 1/2 (q)P] i2) + [Cl/2(q)P 



1(4) 



For the sake of brevity, the derivation of this expression and the coefficients are given in Appendix O!. At this stage 
we have thought it interesting to study qualitatively the influence of the derivatives of (q) in this formalism. For 
this we check the accuracy the approximation used by calculating ^,1/2(9)^.1/2(9) which we write in the form: 



^,1/2(9)^,1/2(9) = 1 + u(q) + R(q)- 

The remainder R(q), whose calculation is given in details in the Appendix OB is in fact an operator 

R(q) ~ A(NW (q)) + l [C(N {2) , A> (2) ' , A (2) ") N {2) (q)P] (2) 



(33) 
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The symbol ~ is used to indicate that symmetric operators [.P]^ n \n > 2 are neglected in the expression of R(q). 
This operator, which should be zero if the moments were constant, provides some quantitative information on the 
accuracy of our approximation as function of the first and second derivatives of N^(q). In Fig. [12] we give a plot of 
the coefficients A and C over a wide range of deformations between 30 b and 350 b. It is seen that sup|^4| = 1.9.10~ 2 

and sup|C| = 9.8. 10 -2 over the whole range of deformations considered here, which suggests that this ansatz is a 

reasonable zeroth-order approximation, despite the large variations in the moments iV^ 2 ). On the other hand, the 
expressions for the coefficients A and C given in Appendix C tell us that A and C vanish independent if one can 
neglect the second derivative of N^ 2 '(q). In that case, J4, 1/2(9) becomes an "exact" solution consistent with the 
truncation we made in the expansion of the kernel overlap. Notice also that such an approximation takes into account 
all the derivatives of N^(q) up to second order and still includes the first derivative of N^ 2 \q). 

0.1F 1 1 ' 1 ' 1 ' 1 



• A 

0.05 - A C 



0- »•••••«•• • 1 • • 

* A A 
▲ 

▲ 

-0.05 - 



-0.1 



100 200 300 400 

q(b) 

Figure 12: A and C coefficients in Eq. ([55)) as function of the deformation. 



The same discussion applies to the operator </a, -1/2(9) which is written in the form 

4-1/2(9) = A_ 1/2 (q) + [B_ 1/2 (q)P] (2) + [C_ 1/2 (q)P] (4) , 

the coefficients A_ 1 / 2 (q), B_ 1 / 2 (q),C_i/ 2 (q) being given in Appendix [C^. In particular, if we neg lect NW{q) this 

operator satisfies Ja, -1/2(9) Ja,+i/2{q) — 1 an d more importantly it is also solution of Eq. (|3ip up to second order. 
Inspection of the equations reveals that the coefficient C-i/ 2 (q) is not completely defined at this order. A complete 
determination of this coefficient is achieved by canceling the contribution of fourth order in Eq. f|3 1[) . Its value is then 
given by 

c-i/aG) = |(^ (2) (9)) 2 [i - ^««) + ^em 

with £(g) = (A>( 2 )'(g)) 2 /A>( 2 )(g). 

To obtain this expression, the third and fourth derivatives of CL 1/2(9) are set to zero. It turns out that these conditions 

are satisfied by the solution given above since we set to zero the derivatives (N ( ^(q)) ,p>2. Thus we conclude 
that it is a solution in the whole range of deformations consistent with our assumptions on the derivatives of N^ 2 '(q). 
To summarize, the operator 

J-1/2® = l-J[^ (2) (9-)(l-^(^ (2) '(9-)) 2 )P] (2) 

+ |[(^v (2) (9)) 2 (i - f m + ^ 2 (9))p] (4) 

is solution of Eq. ([3Tj) up to forth order in the symmetric operator [.P]^- 1 . This expansion takes into account not 
only the derivatives of N^(q) up to second order but also the first derivative of N^(q). Even though it is not 
complete, since we do not include the second derivative of N^ 2 \q), this operator represents a significant improvement 
over what was generally used in previous works. We now proceed to the iV-dimensional case by following step by step 
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the procedure described in this section. Consequently, we will also neglect the contribution of the second derivatives 
of {q) in the expression of the inverse overlap kernel. One finds a justification for this approximation in the fact 
that the general structure of the Schrodinger equation we derive does not depend on the precise determination of the 
coefficients entering in the definition of our ansatz for the inverse. 



B. Determination of the inverse in the iV-dimensional case 



A similar analysis to the one-dimensional case can be performed here by considering instead the operator 

m = 1 (r/v(%)pl (1) + lr/y(%) P i( 2 >) 1 

First it is written in terms of normalized moments according to the expression 

u(q) = [W N (q)P] {1) + ^[NW(q)pf ) +a (0) (q), 
where we have introduced the quantities: 

W N (q) = iV (1) (g)+iC (1) (jV (2) ,jV (0) ) 

The expression of (N^\ N^) is given in Appendix [ClL . According to the definition of J(q) (Eq. (jTDJ) ) this 
operator becomes 

J(q) =I + a (0) (q) + [W N (q)P] (1) + \ [N™ (q)P\ (2) . 

The difference with the one-dimensional case is the occurrence of an additional first-order symmetric operator and the 
fact that all operators are matrices. At this stage we proceed as we did in the one-dimensional case by introducing 
renormalized zero-order moments and using an iterative procedure. Let us mention here that our numerical result tell 
us that the matrix elements of atQ\(q) never exceed 2 x 10 _1 . Consequently a series expansion with few terms gives 
the operator (I + a^(q))^ 1 ^ 2 with excellent accuracy. We recall that we keep the same notation N^°\q) to denote 
this renormalized moment. With these renormalizations we are finally led to solve Eq. (|31|) which can be written as 

J-i/ 2 (q)[I + u(q)]J-i/ 2 {q) = I 

u(q) = [W N (q)P] W + \ \m (q)P] (2) . (34) 

One could attempt to solve this equation with an appropriate ansatz which, in the N-dimensional case, should also 
include odd-order symmetric operators with an ansatz of the form 

J-l/M = A_ 1/2 (q) + [D_ 1/2 {q)P] {1) + [B- 1/2 (q)P} {2) 

+ [E-l/2(q)Pf ] + [C- 1/2 (q)P] i4) . (35) 

This expression should be inserted in Eq. (|34[) but for reasons given in the one-dimensional case, it is advantageous 
to start with an approximation to Ja, -1/2(9) chosen as the series expansion of (/ + u(q))^ 1 ^ 2 up to second order in 
iiiq), i.e. we set 

1 3 

^4,-1/2(9) = 1 - ^u{q) + g« (q), 

with the operator u(q) defined this time by equation (|34[) . It is straightforward to express J_i/ 2 (q) in a form similar 
to Eq. ([35f and to rewrite Eq. ([35f as: 

u,m = i-\[wmp] {l) +-\[^ 

- ^(N {2Y (q)) 2 -l(W N (q)) 2 



3i 

[W N [q)J\ 

I (3) , /^nlW 



16 -(^(g)iV( 2 )'(g-))-)P] (2) 



+ [E_ 1/2 (q)P] w + [C_ 1/2 (q)P] { 
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with the notation 

(X,Y)± = XY±YX. 

The first derivative of Wn(q) has been neglected since it is generally very small. The operator defined by the first 
three lines of this expression is solution to Eq. ([34]) up to second order. The coefficients £^-1/2(9) and C_i/2(g) can 
then be determined by canceling the third and fourth order contribution in Eq. (|34|) . For practical reasons we neglect 
their contributions and derive in the next section a Schrodinger equation with the inverse defined by 

- ^(W N (q)N {2y (Q))-)P] {2) - (36) 

Otherwise the expression would have contained an exceedingly large number of terms. Again we justify this approxi- 
mation by the fact that it does not affect the general structure of the Schrodinger equation that we propose. 

V. DERIVATION OF A SCHRODINGER-LIKE EQUATION 

According to Eq. (fT9f and taking into account the fact that the inverse operator J_i/2(<?) given by Eq. (|36[) is 
Hcrmitian, the contribution of the Hamiltonian kernel to the Hill- Wheeler equation reduces to the form 

i - 1,2<5) 7F!Sl (fl( ° l( ' ) + [,^<1,(, " )Pl<,, 

Now, we proceed as we did in the case of the overlap kernel by introducing this time normalized moments of the 
Hamiltonian defined by 

gto(g) = — -L— ff(")(g)- ' 



and rewrite the equation given above as function of these moments 

^-i /a («)(^ (0) (9) + ^ 1) (^ (1) .^ (0) ) 

+ \[H^{q)Pf ] )J_ 1/2 {q). (37) 

The next step is to extract a Schrodinger- like equation following the procedure described in (23| . To this end we 
insert the definition from Eq. (|36[) in Eq. (|37[) and make an expansion in terms of symmetric ordered product that 
we truncate at order two. We illustrate in the following how we achieve this goal by rewriting Eq. (|37[) in condensed 
notation 

J-i/2{q)W-{q)j-i/2(q), 

with 



J-i/2(q) = T,U[ki)(Q)P] 



where the quantities jii)(q),h^(q) are obtained by identification of these operators with those defined in Eq. (|36[) 
and Eq. ([3T)l respectively. With these notations we can rewrite Eq. (|3"T|) in the convenient form 

p,q,r— 
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which has the advantage of showing the typical terms occurring in this formalism. The expansion we are looking for 
is then derived very simply by means of general formulas that are given in Appendix [X] (Eqs. (|A3|) and (|A4[) ). This 
derivation is straightforward, but requires nevertheless quite cumbersome manipulations that we do not describe here. 
Thus we limit ourselves to giving here the final expression of a Schrodinger equation which takes into account the 
coupling between collective and intrinsic excitations. This equation is written as 

(H(q)-E)g(q) = 0, 

H (q) = S(q) + [T(q)P] W + [U(q)P] (2) (39) 

where the quantities S, T, and U are matrices whose expression is 



h (0) + g( /l (0)J(2))+ + J^j(2) h U j (2) ~ i^SoJid))- 



T = 



1 1 ■ l(1) 1 / ■ l(2) 

- gU(D>)J(i))s - Y6 (j(2) h Wvv A 

+ l(j(2)h ( ^ ) j (1 ))A- ^0'(2) ft (0)i(l))^ ( 4 °) 

For the sake of convenience, the dependence in q has been removed in all the quantities in this formula where we use 
the compact notation 

(XYZ) S = XYZ + ZYX 
(XYZ) A = XYZ-ZYX. 

Furthermore, the expression for the matrix elements of the moments of H(q) suggest that we make the same 
assumptions as for the derivatives of the overlap moments. Thus in the derivation of the formulas proposed here it 
has been assumed that: 

'h$M = for p>3 

hB(q) = for p>l (41) 
h¥jq) = for p>2 

These assumptions serve to reduce significantly the number of terms in the expressions of the different contributions 
S(q),T(q),U(q). 

We discuss now some aspects related to the properties and the general structure of the Schrodinger equation defined 
by Eq. (f3"9")l and Eq. (|4T)|) which contains all the necessary ingredients to build it. First, we note that this Hamiltonian 
is Hermitian and time-reversal invariant. This is easily demonstrated using the properties of the matrix elements of 
the moments listed in Section IIII C 31 As a consequence, if we treat collective and intrinsic excitations on the same 
footing, there is no dissipation but only a reversible exchange of energy between these two kinds of excitation. Notice 
in particular that a linear term as it occurs in Eq. (|39l) does not induce a dissipative force as introduced in studies 
of fission using a macroscopic-microscopic approach [19]. However, It may be that for practical reasons we need to 
treat explicitly a reduced set of intrinsic excitations and take into account the remainder implicitly because they 
are too numerous and their spectrum too dense. In that case we may be led to introduce a dissipation mechanism. 
We do not consider for the moment applications at high excitation energies where such situations would occur (34[. 
Concerning the general structure of the Schrodinger Eq. (|39|) . we analyze separately the three different terms which 
are present in its definition. In Eq. (|40[) we have arranged the terms in such a way that even moments come first 
and odd ones at the end in S(q) and U(q). We recall that odd moments occur only in the multi-dimensional case. 
Accordingly in evaluating S(q) and U(q) in the one-dimensional case we must consider only the first three terms in 
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S(q) and the first six in U(q). Their expressions are readily obtained with the formulas given previously but they 
are very long and for the sake of brevity we do not give them here. Notice that if we neglect the derivatives of the 
overlap moments we recover, as we should, the simple Hamiltonian given in reference [23j . In the multi-dimensional 
case all the terms in Eq. ([H))) must be retained. They involve products of matrices of the moments of the overlap and 
Hamiltonian kernels, and their derivatives. Each of them can be interpreted as effective vertices inducing transitions 
between excitations. At this point it is important to note that these vertices, by the definition of the moments, are 
the result of an average over a wide range of deformations around q and consequently the coupling between collective 
and intrinsic degrees of freedom is very non-local in the present work. This contrasts with the other approach (35| 
which couples the nuclear collective dynamics to internal excitations defined at every deformation along the adiabatic 
fission path. To conclude this analysis it is appropriate to recast our Hamiltonian Eq. (|39p in a familiar form. With 
more conventional notations, first we express its diagonal matrix elements which leads to the expression 

Hu(q) = l[(j^)uPf ) + V ii (q), 

We have taken into account the fact that the matrix T is zero along the diagonal and identified a collective potential 
as: V(q) — S(q). Also we have defined the mass as: Af~ 1 (g) = 2U(q). We recognize standard Hamiltonians which 
describe separately the collective dynamics on each potential energy surfaces built with the various 2-qp excitations 
considered. In the same manner we introduce the coupling between 2-qp surfaces by considering the non-diagonal 
matrix elements of Eq. (|39[) . They are written in the form 



with the definition of a non-diagonal mass as: 




= 2Uij(q). 



This form of the interaction term is interesting because it separates the coupling between collective and intrinsic 
degrees of freedoms into terms of different nature. The first two, which depend on the collective momentum, can be 
interpreted as a dynamical coupling while the last contribution is a potential coupling. In appearance the Hamiltonian 
derived here has a general structure very similar to the classical Hamiltonian given by the authors in Ref. [35| . 
However their derivation and ours rely on very different approaches. In particular the one presented here is based on 
a quantum description from the start, while the classical Hamiltonian in Ref. [35] needs to be quantized after the 
fact. 

For completeness, we could have included in our discussion the explicit calculation of the moments associated with 
neutron- and proton-number operators, in order to use the method (2f| which incorporates the conservation of the 
average particle number in the GCM formalism. Similarly, we could have included the moments of the operator 
associated with the subtraction of the center-of-mass energy of the nucleus. These contributions are straightforward 
to calculate, but are omitted here in order to avoid further complicating the expression for the Hamiltonian above. 



VI. CONCLUSION 



In this paper we have presented a theoretical framework, the Schrodinger Collective Intrinsic Model (SCIM), 
which allows in a microscopic way the simultaneous coupling of single-particle and collective degrees of freedom. 
Such an approach is based on a generalized GCM, where the general GCM ansatz of the nuclear wave function is 
extended by a few excited configurations. In fact, one considers as generating wave functions not only HFB ground- 
state configurations with different values for the collective generator coordinate but also 2-qp excited states. Such an 
approach has the advantage of describing in a completely quantum-mechanical fashion, and without phenomenological 
parameters, the coupling of qp degrees of freedom to the collective motion of the nucleons. Our derivation of the 
SCIM proceeds in two steps. The first step is based on the generalization of the symmetric moment expansion of the 
equations derived in the framework of the GCM including the coupling between collective and intrinsic variables. The 
moments occurring in the equations are matrices whose dimensions depends on the number of excitations introduced 
in the description. Contrary to the usual case - without qp excitations - all moments, odd and even, must be included 
in the summation. Such an exact expansion in terms of local operators of the generalized Hill- Wheeler equation is then 
transformed into a local Schrodinger equation by inverting the expansion of the overlap kernel. Then a second-order 
differential Schrodinger equation is derived by limiting the summation to second order terms under the assumption 
that the series expansion of the overlap and Hamiltonian kernels converge rapidly. The derivation of the inverse of the 
overlap has been seriously complicated by the fact that the moments and their derivatives vary rapidly as function 
of the deformation, as observed in the numerical study of the overlap kernel in a wide range of deformations in 236 U. 
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A quantitative information on the accuracy of our approximation as function of the first and second derivatives of 
second-order moment of the overlap has been given in the present paper for the scalar case. The Schrodinger equation 
is finally written in a convenient form that has the advantage of exhibiting typical terms occurring in the formalism 
such as, for instance, the potential coupling, and the dynamical coupling between intrinsic and collective excitations. 
There are a number of avenues that could be pursued now with this new formalism. Among them, the study of the 
fission dynamics and the coupling between intrinsic and collective excitations in the descent from saddle to scission 
represents one of the most challenging problems in many-body theory. Also, a particularly interesting use for this new 
approach is in those nuclei which exhibit excitation spectra showing collective as well as non-collective effects. This 
would surely improve the energies of the states, which is one of the problems of collective models today. In all these 
studies, the microscopic nuclear wave function will be analyzed to determine the influence of the different degrees of 
freedom on nuclear properties. This would of course add to the computational burden, but the ingredients to perform 
the calculation are ready for the most part. The derivation and calculation of energy kernels will be presented in 
detail in a separate publication. 
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Appendix A: SYMMETRIC ORDERED PRODUCT OF OPERATORS 

The symmetric ordered product of operators (SOPO) is defined by [23[ (p. 421): 

[A(q + s/2,q-s/2)P] {n) = 

n 

— C«P n - q A(q + s/2, q s/2)P« (Al) 

9=0 

where A is any operator and P is the derivative operator P = jJL For instance, the first three SOPO are: 

-{AP + PA) 

-(AP 2 + 2PAP + P 2 A) 

4 

where the action of PA on a vector g(q) is: 

PAg{q) = A^g(q) + APg(q) (A2) 

with AW defined as 4" = (P n A) — (i)™^-^. Let us note that for a Hermitian operator A, A^ is Hermitian or 
anti-Hermitian according to the parity of n: 



[AP] 


(0) 


= A 


[AP] 


(i) 


1 

2 


[AP] 


(2) 


1 

4 



(a^Y = (l)"iW 



For any operators A, B, and C, which depend on the collective variable q, we have the following properties: 
- The symmetric ordered product is linear: 

[AP] (n) + [BP] (n) = [{A + B)P] (n) 
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The symmetric ordered product preserves the hermiticity: 

A+ = A ([AP] {n) y = [AP] (n 

{[AP] (n) [BP] iq) Y = [BP] (q) [AP] in] 



f + Z A B => (\AP] (n) \BP] (a " 



- More generally, the product of two and three SOPO can be expanded as a linear combination of symmetric ordered 
product. This property is used throughout the paper in the derivation of various expressions. 

We give these expansions in a form convenient for retrieving the contribution of a SOPO of a given order 

n+q min(n,i) 
i—0 s-max(0,i-q) 



[ i (i-») B W P ](«+H) (A3) 



[AP]^[BPf\ C P]^ = E ^ E E 

p=0 i=0 t=max(0,p- q— r) 

(n+q+r- 



C* C&ZK-lY [A^-^F(B, C, r, i, t)P] {n+q+r ~ P) (A4) 



with 



(*.«) . (t) 



F{B,C,r,i,t) = E (B (J - S) C (S) ) 

s-max(0,i-r) 

In particular, using (|A4|) with n = 0, g = 1, r = we find: 

id) , 1 



A[BP] {1) C = [ABCP] W + ±(abcM - A^BC 



A[BP] (2) C = [ABCP] (2) + [(ABC™ - A^BC)P] (1) 

- ( A&BC - 2A«BC« + ABC^ 
4 V 



(A5) 



4 

Formulas (|A5[) and (|A4|) are the formulas used to develop the formalism in terms of normalized moments. In fact 
after setting 

1 



we obtain 



A = C-- 

%/aW 

—^\BP] {1) —L= = \BP} {1) + -C^\B,N^) 
1 [BP] [2) -L= = [BP} {2) +iCW(B,NW) 



C-(2)( BjiV (0)) ( A g) 



with the definitions 



C W(A,nW) = —^=A ( -fLJ)' - ( -rL=)'A ' 



c( 2 )(A,iv( )) = i[(_L-Y 



r i y^ i 

I 1 1 



1 1 •"' 
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Appendix B: OVERLAP KERNEL 

The matrix elements of the overlap kernel are expressed as Nij(q,q') = ($i(q)\$j(q')), where is defined by 

Eq. (EU). 

The calculation of the overlap kernel between two HFB states at different deformations is detailed in [2H [27j. We 
use the formalism i n 12a. 12711 b ecause we are most familiar with it, but of course the same formulas exist elsewhere in 
the literature (e.g., |ldL l23l. |25|| and references therein). The generalization to 2-qp excitations is straightforward. 

Let us introduce the quantities "S, T and Y" of H&EzI 



Sij 



T ■ 



(*(9)l*(g / )) ' 

($(q)\at q 'ap'mq')} 

($(q)Mq>)) 
$>(q)\afa<\Hq')) 



^ ~ <<%)|<%')> (B1) 
where the particle operators {a q ,af q } are related to the qp operators {fjf,??^ 9 } through the BCS transformation: 



Vi = u q a q - v q ap 

-\-q q +0 . q q 

\ = u i a i + v i a l 



(B2) 



For k, I > 0, and k, I their time reversed states, we have: 

v(u'>'(6 ij -S ij )+v(vfY j i)], 

<<%)!«<%')) = {Hm{<i l ))Y, T iv' T lt\ u >i Y pp' 

pp' 

- u q v q k (5 pp > - S pp >) + vfulS p ip - vfv q k Tpp> (B3) 
where the transformation r qq is defined through the relation 

E^V" ( B4 ) 



4 q 



The different quantities in Eq. (B3) can be expressed as (2H |27|: 

($(<?)|$(g')) = det(r qq ')det(Z) 

T fj = (u*'z- 1 v*T<">') ji 

Y fj = -{T qq '- 1 u q z'- 1 v q ') l3 (B5) 
with Z = u q {T+)- 1 u q ' + v q rv q ' and Z' = u q ' (t)- 1 ^ + v q 'r + v q . 

Appendix C: Calculation of I + u(q) and related quantities in terms of normalized moments 

1. Calculation of I + u(q) 
According to Eq. (fT6|) . the operator I + u(q) is 

7+ " (5) = ,+ 7mw) (lNm(s)p]m 

+ > m ^vm> (C1) 
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By use of Eqs. (|A5[) and (|A6I) . can be expressed, after straightforward calculations, as 

u(q) = C+(q) + [W N (q)P] {1) + l[N( 2 Hq)P] i2) 



with 



and 



W N {q) = N { - 1) {q)+iC {1) {N i - 2 \N^) 

C+(q) = t C W(N( 1 \NW)-~cW(NW,NW) 



(C2) 



C«(A^ (0) ) 



=A 



A- 



\/N^Hq) J « 



C( 2 )(A,iV( )) 



1 



A- 



1 



V?m(qjJ ^/NW(q) 
1 \ . / 1 



.4 



It is worth noticing that in the scalar case, the coefficient C^(A, N^) and N^(q) are zero so that the operator 
u reduces to 



l(q) = -~CW(N<V,N<V) + ^[NW(q)P] {2) 



2. Calculation of the operators (I + u(q)) 1/2 and (7 + u(q))~ 1/2 

For sake of simplicity, we restrict our discussion to the standard -one dimensional- case. Ji/ 2 (q) the square root of 
the operator 1 + u(q) is approximated by the first terms of the series expansion up to the order u 2 (q): 



Jx/M = {! + u{q)) 1/2 = 1 + \m - \u\q) 



(C3) 



with 



m = \[N^(q)Pf } 



Using the property (A.l), (u(q)) 2 is 



u\q) 



^{N^\q))\l[(N^\q)N^(q) 
2(N^' (q)) 2 )p] i2) + [(N^(q)) 2 P] i4) 



and the operator Ji/ 2 (q) can be expressed in terms of SOPO as: 

Jl/2(q) = A 1/2 (q) + [B 1/2 (q)pf } + [C 1/2 (q)P] (4) 



(C4) 



2G 



with 



A 1/2 (q) = l-^(NW(q)) 
The same derivation applies for J-i/2(q) the inverse of the root square of 1 + u(q), which can be expressed as: 



with 



j-i /2 (q) = i - \m + lu 2 (q) « (/ + mr 1/2 

(2) 



A-y 2 (q) + [B_ 1/2 (q)P] ' + [C-i/ 2 (q)P] 



,(4) 



A_ 1/2 (q) = l-±(NW"(q)) 

S_i /2 (g) = 

[c_ 1/2 ( 9 -) = i(iv( 2 )( g -)) 



(C6) 



3. Validity of the approximation used to calculate the operator (I + u(q)) 



1/2. 



The validity of the truncation made on the expansion of the square root of the overlap kernel Ji/ 2 (q)i Eq. (|C3|) . 
is checked here, by estimating the extent to which the relation (/ + u(g)) 1//2 (7 + ii(q)) 1 ^ 2 = (I + u(q)) is satisfied. In 
other words, using the expression of (I + ii(q)) 1 / 2 given in Eq. (|C3[) . the residual R(q) is defined by the relation 



(i + u(q)) 1/2 (i + m) 1/2 = (i+ m) + r(q) 



with 



R(q) = -l(u(q)) 3 + ^(u(q)r 

Consistently with the truncations made in section [Til the residual R(q) is restricted to its lower orders in its develop- 
ment in symmetric ordered products. Thus R(q) is 

R{q)^A{q) + \[C{q)N^{q)P] {2) 

Neglecting the derivatives (N^ 2 \q))^ for p > 2, after tedious but straightforward calculations using Eqs. (|A3|) and 
14[) , the coefficients A(q) and C(q) are 



A(q) = ^(N^"(q)) + ^(^ 2 )"( g -)) 
C(q) 



2 ; » 



1 ( si m'{q)(m{q)X 

" 9 (2) ®) +32 — wm J 



\ 



5 ( (^ (2) "(9)) 2 (^ (2) '(9)) 2 (2Y , v8 

+ — 34-^ i— ^ '— + ( iV (2) (qY 

213 I JVW(g) V W < 

The smaller the residual R(q), and by extension A(q) and C(q), the best the approximation. In particular, we find 
A(q) = C(q) = for N^" (q) = 0. Similar conclusions are found for J -1/2(0) • 
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